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Abstract 

A simple slave-boson representation combined with the Hartree-Fock approx- 
imation for the Hund's rule coupling is introduced for a doubly degenerate 
narrow band, which bears a direct relation to that introduced previously in 
the nondegenerate case. Namely, one keeps the fermion representation of the 
spin operator to recover properly the energy of fermionic quasiparticles in the 
presence of an applied magnetic field. A simple two-parameter mean-field 
analysis of the metamagnetism is provided, with the emphasis on the role of 
the Hund's rule coupling. We also analyse the appearance of the spin-split 
effective masses in the applied field and for nonhalf-filled-band situation. The 
Mott-Hubbard boundary is determined at nonzero temperature (T > 0); it 
shifts towards lower interactions with increasing T and the field signalling 
the precursory localization effects, explicitly exhibited in the behavior of the 
magnetic susceptibility calculated in the Appendix. We also formulate a more 
general two-parameter rotationally invariant approach for an arbitrary degen- 
eracy d of equivalent orbitals and show that the Mott-Hubbard transition at 
zero temperature and at any integer filling n > 1 is always discontinuous. A 
brief overview of experimental situation is also made. 

I. INTRODUCTION 

In recent years one observes a renewed interest in the metal-insulator transition of the 
Mott-Hubbard type, with a particular emphasis on the effect of orbital degeneracy and 
associated with it Hund's rule coupling, both appearing as additional factors This 
interest is induced by the new works on perovskites 0, NiS2-a;Se x 0, as well as on the 
canonical system V2-^03 [|. In most of these compounds one observes an anomalous 
metallicity even in the paramagnetic state Among them is the metamagnetism and 

a relatively strong applied magnetic field dependence of thermodynamic properties ||, in 
addition to the field induced metal-insulator transition || . These phenomena of localization 
occur in compounds of various crystal structure and even in the same magnetic phase, as 
in the case of NiS2-o;Se x . Therefore, it seems that the underlying microscopic mechanism 
is rather general, neither strongly dependent on a particular crystallographic structure, 
nor on the form of the density of states near the Fermi energy and the type of magnetism. 
Additionally, electronic properties of the system NiS 2 - 3; Se :r are well described by considering 



a half- filled e g band @ composed of 3d states (due to Ni 2+ ions), which hybridize with 
completely filled p states due to S 2_ or Se 2 ~ (the last states are believed to play only a 
passive role near MIT). 

The purpose of this paper is to carry out a simplified analysis for a system composed of 
orbitally degenerate but otherwise equivalent band states that nonetheless addresses some of 
the principal phenomena in an applied magnetic field and at nonzero temperature. Thus, one 
can examine the system behavior as a function of experimentally controllable parameters. 

The analysis of Mott-Hubbard localization in the nondegenerate band case was based on 
the Hubbard model and provides a continuous transition at T = and a discontinuous 
transition at T > ||. These results were subsequently confirmed in the limit of infinite 



dimensions ||10|| . This means that this quantum phase transition possesses an upper crit- 
ical dimension and the Gutzwiller approach represents a correct mean-field theory. Also, 
close to the localization, depending on the band filling, one observes metamagnetism || or 
metamagnetic behavior J7|,rT|, spin-split effective masses [|TT| , |T2"| , and a transition from an 
antiferromagnetic semimetal to an antiferromagnetic insulator The related Gutzwiller 
approximation scheme for a doubly degenerate Hubbard model provides a discontinuous 
transition already at T = a change induced by the Hund's rule coupling. Obviously, 
the Hund's rule should make the localization easier (i.e. diminish the critical value U c of the 
intraatomic interaction), since it favors energetically the high-spin atomic state against any 
normal Fermi-liquid state. Hence, we can separate the effects associated with the Coulomb 
interaction from those due to the Hund's rule coupling. 

The structure of this paper is as follows. In the next section we introduce a slave- 
boson approach combined with the Hartree-Fock approximation for the Hund's-rule term, 
which reproduces in a simple manner the similar results of the full slave-boson analysis 
by Hasegawa [[I] in the paramagnetic state for zero field and temperature. Additionally, 
this representation provides a correct expression of the energy a fermionic quasiparticles in 
the applied magnetic field. Within this scheme we determine the magnetic susceptibility, 
metamagnetic properties, spin-dependent effective masses, as well as determine the Mott- 
Hubbard boundary at nonzero temperature, and its shift towards lower temperatures in the 
applied field. The similarities and differences with the nondegenerate-band case are stressed. 
In the second part we propose a simple rotationally invariant formulation of the problem 
for an arbitrary orbital degeneracy of electron states and compare it with the Gutzwiller 
approach. Within this approach, a difference between the Hund's rule and intersite magnetic 
coupling appears naturally, although the last point is not discussed in detail in the present 
paper. The correlated metal-paramagnetic insulator transition takes place at any integer 
filling n and is discontinuous for n > 1. 
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II. THE PARA- AND META-MAGNETIC STATES OF DOUBLY DEGENERATE 

HUBBARD MODEL. 



A. The slave-boson method combined with the Hartree-Fock approximation for the 

exchange term. 

Doubly degenerate Hubbard model in the situation with equivalent orbitals is described 
by Hamiltonian: 

ijla il ^ ' i i il 

(1) 

where au a is the annihilation operator of electron from lattice site i on orbital /(= 1, 2) and 
with spin cr(=t, j), n %ia is the particle number and Su = (S£, S~ n S* t ) is the spin operator, tij 
is the hopping integral, assumed for simplicity the same for the two orbitals (we choose also 
tu = 0). Parameters U, U and J, are respectively, the intraband Coulomb, the interband 
Coulomb and the Hund's-rule-exchange integrals. For e g band we have: U = U — 2J. 
Finally, in the last term h = ^gfi B H a , where H a is the applied magnetic field. Because of 
the vanishing orbital moment in the first order, the magnetic field affects the spin degrees 
only (we neglect also the Landau quantization effects for low fields). 

We make the mean-field approximation for the exchange term. Likewise, we decouple in 
the same manner the interorbital term, since U — |J = U — |J. Therefore, we have 

nnn i2 = nn(n i2 ) + (nn)n i2 - (nji)(n i2 ) , 
Sfi • S i2 = Sf^Sfz) + (Sf^Sa - (5 , i z 1 )(S , ^ 2 ) , (2) 

where (A) = Tr(Ae~ l3HH - F ) is the average value of operator A with Hamiltonian TCh-f'- 

ijlcF il il ila ilcr 

— -KNn 2 + - JNm 2 - /iN e . (3) 
4 8 

In this equation we have added the chemical potential part (—fiN e ), with N e being the 
number of electrons in the system. The n and m are defined by equations: 

(n u ) = -n , (4) 
2<fiS> = \m , (5) 

and is the number of sites, K — U — \ J. The mean values of tin and Sfi are independent 
of site and orbital indexes because of translational symmetry and the equivalence of the 
orbitals for half- or nearly-half-filled band configurations (in this part we consider only the 
limit of n equal to or close to two). 



We employ the slave boson representation of Kotliar and Ruckenstein [13| in the form: 
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IM> = 4l w > > 



(6a) 
(6b) 
(6c) 



where we introduced new boson-fermion vacuum \v) instead of the Fock-space vacuum 
{|0, il)}. So, we have added only the orbital index to the original formulation. The op- 
erators el, pt and d\ are boson operators and ft is fermion operator. This representation 
must be supplemented by the constraints 



and 



Qila = fuJila ~ PilaPilc ~ d\,du - 



P a = e^eu + J2p] la pua + d^da -1 = 



(7) 



(8) 



Relation between original (a« CT ) and new (fua) fermion operators is: 



(9) 



The constraints given by eqs. (0) and (j8|) can be enforced by Lagrange multipliers 



(i) 



(2) 

and \ u ^. We obtain following Hamiltonian: 

7~L t ot = 7~Lh-F + ^con 

1 



where 



and 



^ 'tijzllaZjltrfilafjla + U ^ ^il + 
ijla il 



-Kn — n) — a(h + —Jm) 



filer file + ticon + C 



c 



■\KNn 2 + - JiVm 2 , 
4 8 



(10) 



(11) 



(12) 



Note that using the fermion representation for the spin operator = | (fj^fu^ — fjnfiii), 
we obtain correctly the Zeeman term and also, the Hund's rule coupling in the form 
of a molecular field. Equivalently, one could utilize the Bose representation Sf t = 

\ (puiPui ~ PuiPiiij ) but then only the z-component is well defined and the result in the 
saddle-point approximation does not provide the quasiparticle energies with the proper Zee- 
man term (see the discussions of this point in Ref.6b, where the full spin rotation invariant 
form of the SB representation is invoked). 

The factor zu a is not unique, but when we make the choice JT"3||: 



%ila 



(13) 
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we obtain correct result in the saddle-point solution. 
The partition function is given by: 



Z= D[e,J}D\p,p^D[d,S]D[fji}]Jd\VdX$lex P - C(r)dr 

J ila L ^° 

where (3 = 1/ksT is the inverse temperature, 

C ( T ) = ^2fL (d T -]l~a( y h+^Jm) + Ag) 6 tj + UjZ^Zjy 
ijitr L v J 

+ E 4 (* + \? ; ) e, + E?l (a + 4" - a») Kfa + 

il ila 

+ E 4 (a + a«;> - £ a£ + c/) - y, ^ + c . 

il \ a J il 

is the system Lagrangian, and 



(14) 



fjla + 



(15) 



A 4 — M — 



(16) 



is the effective chemical potential. 

Next, we make the saddle-point approximation, i.e. assume, that Bose fields and La- 
grange multipliers do not depend on time r, number of site i and orbital index I, so e^(r) — > e, 
Piia( r ) Pa, du(r) — > cf, A^J — * Act 2 "* and Ay — > A. To obtain saddle-point solution we must 
minimize free energy with respect boson fields and Lagrange multipliers. From partition 
function (14) we obtain following expression for the free energy: 



F = -k B T^\n 



Ma 



1 + exp 



M - E] 
k B T 



1 



-NJm 2 + 2NUd 2 + 



-2N*£\W (p 2 + d 2 ) + 2iVA« ^ e 2 + ^> 2 + d 2 -lj + 



— -NKn 2 + uNn , 
4 

Fermi quasiparticle energy is given by: 



E ka = q a e k - a{h + - Jm) + A^ 2) , 



where ek is bare band energy, q a = z 2 and 



(ep a + dpaf 



(l-d2_ p 2 )(l _ e 2_ p |) 



(17) 



(18) 



(19) 



Thus, one encounters here two molecular fields: one coming from the Hund's rule and 
the other (a| 2 } - Xf coming; from the electronic correlations. 



Differentiation of the free energy with respect Lagrange multipliers provides the mean 
field version of the constraints: 



<>l pl + d 2 = ±-J2n ka = \n a , (20) 



m ~ u ^ ^ u ~ n ^ kCT ~ 2 

dF 



dXt k 



where 



=► i- e 2 -E^- rf2 = ' ( 21 ) 



^ l + exp|)9(£W-AZ)] (22) 

is Fermi-Dirac function. From eqs. (§) and (0) we obtain that (n^) = \n a — (n + am)/ A. 
Instead of computing the derivatives with respect to the remaining fields we, put e and p a 
expressed via d, n and m into free energy ([[7]). Next, we introduce new variables (3 and /? , 

(2) 

such that = Po + crfi. We obtain the following free energy 



- = -2k B T- Y In [l + exp ( l£ll3A 
N B N ^ [ \ k B T ) 



Quasiparticle energy is given by: 



I i 
+ -Jm 2 + 2Ud 2 - (3m+ -Kn 2 + Jin . (23) 
8 4 



E ka = g CT e k - cr(h + ^Jm - p) , (24) 



and effective chemical potential by: 



/z = /i-A) = /i - -Kn - Po . (25) 



From ([19]) we obtain q a in a form 

( \/l — |n + d?\Jn + am — Ad? + d\/ n — am — Ad? 

1° = 4 Ta \T— S ( 26 ) 

(4 — n — am) [n + am) 

This factor is the mass enhancement and in general case (n ^ 2) leads to the spin- 
dependent effective masses. 

One should justify the approximations made. First, for e g band J ~ 0.2U , and therefore 
U = U - 2.5J ~ 0.5U. Hence one can say that making the Hartree-Fock approximation 
for the interorbital-interaction part is qualitatively correct even for U ~ W. Below we make 
this argument more quantitative by comparing the ground-state energy obtained within 
different approximation schemes. Mixing the slave-boson scheme with the Hartree-Fock 
approximation for the Hund's rule exchange has this advantage that we can single out what 
feature is specifically due to the Hund's rule coupling, as discussed next. 
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B. Case n = 2,T = 0, and h = 0. Effect of exchange interaction on the metal-insulator 

transition (MIT). 



We consider first the simplest situation: a paramagnetic state in the half-filled band case 
and at zero temperature and in zero applied magnetic field 

At temperature T = the distribution function (p2|) takes the form 

n k = e(p-B k ) = {«^^>?' (27) 
1 1 gdy E k ^ fi , 

where = $o e k and $ = Qa — Qw = 8cP (1 — 2d 2 ) . We obtain following expression for the 
ground state energy: 

-W/2 



/w/z 
dep (e) (/I - $ e) 0(71- $ e) + 2Ud 2 + 2p 
-W/2 

= 2%e + 2Ud 2 + K, (28) 



where: 



rW/2 

e = 2 dep(e)en(e) (29) 

J -W/2 

is the mean bare-band energy, n (e) = (/Z — $oe) , and p (e) is density of states. Minimizing 
the functional (p8| ) with respect <i 2 we obtain the solution: 

which coincides with the results in nondegenerate-band case §|[|, i.e. the metallic state 
would become unstable for U = Uq = 8 |e| if the system evolves continuously with growing 
U. However, we shall see, that in the degenerate system this instability point Uq is replaced 
by a critical point U c < Uq, at which a discontinuous transition takes place. 
Ground state energy of the metallic state can be expressed explicitly as: 

Eq 1 = -2|e| - + -U- - J. (31) 
11 32 |e] 2 2 v ; 

In the Mott-insulating spin-disordered state, i.e. for d 2 = 0, we have 

Eq = K — ^-J=U — 3J. (32) 

Equating the energies (|32| ) and ( p8]) we determine the critical value of J = J c and U = U c , 
for which the transition from a correlated metal to a paramagnetic Mott insulator takes place 

J c = 4|e| (l - |0 * (33) 

For J < J c the metallic state is stable. Also, for J > the transition is of the first order, 
whereas for J = it is continuous and corresponds to that for nondegenerate case. For 
example, for J/U = 0.1 we obtain the critical value U = U c = 5.1|e| . The diminution of U c 
below Uq is obviously caused by the Hund's rule favoring the high-spin state on each atom. 
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C. Case n = 2, T = 0, and h ^ : Metamagnetic transition. 

We discuss next the homogeneous magnetized state starting from the paramagnetic state 
at T = 0. From eq. fl26|) we obtain band narrowing factor in this case in the form: 



16rf 

$ = Qo 



4 — m 2 



1 



l-2rf 2 + --y/4(l-2rf 2 ) -m 2 
2 



(34) 



Thus, the effective mass will not depend on spin for the half-filled band case. 
The average occupancy of each orbital state with spin a is 

2 fW/2 i-(jZ+a(h-/3+Jm/4))/ q< 



n °=N 

k 



V /it,, = 2 dtp (e) n CT (e) = 2 / dep (e) , (35) 

W/2 J-W/2 



where n a (e) = 6 (p + a (/i — (3 + Jm/A) — q a e) . We choose the featureless (rectangular) 
density of states: 

(£)= f^for £ e[-iy/2,W2], (36) 
I U otherwise. 

The overall futures should be independent of the detailed shape of p (e) , as we determine 
global quantities (m, Eg, d 2 ) involving integrals over the filled part of the relevant bands. 
Putting above DOS into (|35D we obtain: 



2 r(JC+a(h-p+Jm/£))/® 

n„ = — / de 



W J_ w/2 W$ 



i a 1 t ^ W ® 
h — p + -Jm cr + p + 



4 / r 2 



(37) 



Summation over a of the above equation leads to the condition p = 0, but multiplying 
it by a and then summing over a, leads to the expression for the molecular field of the form 

1 1 

8 = — $Wm + h + -Jm . (38) 
4 4 

Because of $ is a nonlinear in magnetic moment m, the molecular field (3 depends also on 
m in the same manner. Inserting f3 into (p3|) for T = 0, we obtain: 

E M = — — = - (4 - m 2 ) $/8 + 2wd 2 - Am - Ijm 2 + k, (39) 
1\W 8 

where the reduced parameters are: 

u = U/W, j = J/W, k = K/W, h -> /i/W. (40) 

We determine the saddle-point solution minimizing Em with respect d 2 and m. Differ- 
entiation with respect to m leads to 

d 4 (m 2 - 16a) + a (l6d 2 + m 2 - 4) = . (41) 
where a = (h + jm/4) 2 . The differentiation with respect d 2 yields 
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i (u + 4c/ 2 - 1) ^4(l-2c/ 2 ) 2 -m 2 - (1 - 2c/ 2 ) 2 + ^m 2 + 2c/ 2 (l - 2c/ 2 ) = . (42) 

For h = and ii ^ 2 from (|30|) we obtain d 2 = (1 — u/2)/4. The change of direction of 
magnetic field h does not change d 2 ; thus (<9cZ 2 /dh)h=o = 0. Next, differentiating with respect 
h eq. pi| ) and taking into account the above relations we obtain the following expression 
for the magnetic susceptibility: 



. , dm 



4X0 = , (43) 

h=0 * 2+u J 



where Xo = 2/i 2 3 /H / is the Pauli susceptibility (cf. Appendix A). The magnetic susceptibility 
was derived for the metallic state, so from taking into account the condition (|3^) we obtain 
the metallic state for j < (1 — u/2) 2 ^ |~^ 2 ^ 1. From, the last inequality we obtain the 
condition for the magnetization in small fields: 



dm 
m « — — 
oh 



h = -n—, > z r • (44) 

h=0 2+u J J 



This condition reduces solutions of eq. (f4l|) to: 



2 —8a + m^Ja (16a + 4 — m 2 ) 
m z — 16a 



Substituting d into eq. (|42j) we obtain following equation for magnetization: 

-64a 3 + 12a 2 m 2 - 3am 4 /4 + m 6 /64 + 64a 2 w - m 4 w/4 + 
-16aw 2 + 32a 2 w 2 + m 2 u 2 + 12am 2 w 2 + m 4 w 2 /8 + 

-16aw 3 - m 2 u 3 - 4aw 4 + mV/4 = 0. (46) 

The pair of solutions (m,d 2 ) determines the ground state energy (p9|). Ground state 
energy of a Mott-Hubbard insulator in the spin disordered state is: 

— —2h - j/2 + k. (47) 

From the preceding subsection we know that the metal-insulator transition take place for 
j = (1 — u/2) 2 . From the expression (H3T) for the magnetic susceptibility we see that except 



for the case j = and u — 2, there is no singularity in x &t the transition. This means that 
the transition for j > is of the first-order. For j = we have second order transition, as in 
the nondegenerate band case, where \ is singular (cf. Appendix A). 

Field dependences of d 2 and m is displayed on Figs. lab, respectively. Close to transition 
point the system exhibits a metamagnetic behavior i.e. the m (h) dependence in curved 
upwards. With the increasing field h the doubly occupancy diminishes because of the growing 
spin polarization and at the transition point d 2 jumps to zero, while the magnetic moment 
jumps to its saturation value. The magnitude of magnetization jump Am and the critical 
field h c depend on the values of j and u, as shown in Fig. 2. One observes here an example 
of a magnetic-field-induced localization, with a formation of a ferromagnetic (field-induced) 
insulator. The transition from the metal to the insulator can be either first or second kind, 
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depending on values u and j. Regimes of u and j, where the particular kind of a transition 
takes place, are displayed in Fig. 3. We mark the area of insulating phase in zero field bounded 
by the curve determined from Eq. (|33|). Insulator polarizes at T = to a saturated state in 
an infinitesimal field; jump of the magnetization is then equal to Am = 2. Upon increasing 
of j the critical field h c decreases; likewise the magnetization jump. From Fig.2 and Fig. 3 
one could seen that for u > 0.5 a second-order transition does not occur. Obviously, our 
analysis is valid quantitatively only for j substantially smaller than u. Also, the critical field 
for metamagnetic transition becomes experimentally accessible only close to MIT (u — > 2) 
and reduces gradually to zero with growing j. This prediction could be tested experimentally 
(see below). 



D. Case n < 2 and T = 0: spin-split masses. 

Consider now briefly the situation n < 2 when the metallic state is stable. We examine 
the effect of magnetic field on the transition to the saturated state, as a function of the 
parameters u and j. Ground state energy for a polarized state is provided by the expression: 

E M (n) = = — — q a (A — n — am) (n + am) + 2nd 2 — -jm 2 — hm + -kn 2 . (48) 
iV W 16 8 4 

u 

Energy of the saturated state is obtained by taking into account eq. ( p6|) and making 
substitution d 2 — > and m — > n in (1481). We have then: 



Ei (n) = —-n (2 - n) - ^jn 2 - hn+ ^kn 2 . (49) 

In Fig. 4 we display magnetization m, doubly occupancy d 2 and the mass enhancement 
q^ 1 as a function of applied magnetic field for u = 1.98 and j = 0.05, for the band fillings 
n = 1.86 and n = 1.8 . The important feature is that the mass enhancement is spin 
dependent (see the lowest panel). Also, there is a possibility of first and second order 
transitions to a saturated state. Regimes of u and j , where this two types of transition 
are possible for filling n = 1.8 are shown in Fig. 5. By comparing this Figure with Fig. 3 we 
see that regime of stability of an insulator for n = 2 corresponds here to the area lb with 
excluded doubly occupancy (EDO), when the system is in a magnetically saturated state 
in a vanishing field. In the case of one-band model, i.e. j = 0, EDO state is reached for 
u — > +oo in h = 0. For j ^ 0, the doubly degenerate case, this state is achieved for finite u. 

One should note that the effect of spin-split masses arising in the magnetically polarized 
state is not associated with the emergence of spin-dependent density of states in the bare 
band. This can be seen from the fact that the quasiparticle energies E^ a (cf. Eq.(24)) lead 
to the following density of states: 

<1<T \ J 

where (e) is the density of bare states per spin. Thus, the enhancement due to the 
correlations is distinct from the change of the density of states caused by the presence of 
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the effective field. The many-body nature of the enhancement factor can be also seen by 
writing the quasiparticle energy in the form 



E ka = £k- ah + (q a - 1) £ k - a \ jJ m ~ p \ = £\l - crh - ^ {u = £ k ) , (51) 
where the selfenergy is 

E a (cu) = -(1 -q a )uj-a Qjm-/3J . (52) 

Taking into account the well-known definition of the mass enhancement in the Fermi 
liquid 

^=li m (l + ARe£ CT (u;)Y (53) 

we have that m /m* = q a , i.e. is indeed spin dependent. This spin dependence, which 
will lead to the strong field dependence of the linear specific heat close to the metamagnetic 



transition (see Spalek at al. [11]) is characteristic of the almost localized fermions and should 



be determined experimentally. 

III. MOTT-HUBBARD BOUNDARY AT NONZERO TEMPERATURE 

A. The phase diagram: h ^ 0. 

In the preceding Section we considered the system properties at zero temperature. We 
concentrate now on a more realistic case of non-zero temperature in the half-filled-band case. 
The low temperature (Sommerfeld) expansion || of the free energy (|23|) for the metallic 
phase leads to the following expression for the constant DOS: 

= _ (4 - m 2 ) $/8 + 2nd 2 - hm - -jm 2 + k - ^r 2 , (54) 
N W 8 3^? 

where r = fcgT ' jW. No higher-order term in r appears for this DOS. 

Hamiltonian describing the insulating phase can be rewritten in the form 



i i il 

= [Knana - J (S 2 - S 2 X - S 2 2 ) - 2hSf] . (55) 

i 

where Si = Sn + Sj2 is the total spin per site. The eigenstates of Hi are singlet and triplet 
configurations. The partition function is then given by: 



Z = Tr (e-^) = e~ pKN [e^' 2 + e^ 2 + e^ 2 (e 2 ^ + e' 2ph )] = 

= e- pKN [e-^ J ' 2 + e^ 2 (1 + 2 cosh (2/3h))] " . (56) 
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This expression leads to the free energy for the Mott insulating state in the following 
form 



F > 1 i n [ Z ] = k - -L In [e" 3/3J / 2 + e^ 2 (1 + 2 cosh (2(3h))] 



NW NW(3 1 1 Wp 

= k - rln [e" 3j/2T + e i/2r (1 + 2 cosh (2/i/r))] . (57) 

The result in the limit r — > reduces to the expression (0). The magnetization in this 
state is 



m 



2Tr (Sfe-Wi) 



(3KN p pj/2 ^ e 2(3h _ e -20h\ 



e _p KN [ e - 3 /W/2 + e /3J/2 (J + 2 cosh (2/3/i))] 



4sinh 2/3/i , . , 2 

= —^5 / ; ^ = tanh (h r) ^ , 58 

e" 2 ^ + 1 + 2 cosh (2/3h) K/ ' e-^-i + 1 ' v ; 

For j = we obtain m = 2 tanh (h/r), the expression for the magnetization of non- 
interacting spins. Such a situation arises because the saddle-point approximation in its 
essence is a single-site approximation (the intersite interaction arises from the quantum 
Gaussian fluctuations around it). 

The field dependences of the magnetization and of the doubly occupancy are shown in 
Fig. 6. We see that upon increasing j while keeping fixed u the critical metamagnetic field 
is reduced. It indicates, as before for T = 0, a significant role of exchange interaction and 
associated with it Hund's rule. Also, the magnetization curve is slightly curved upwards 
in small fields. The upper part of the m(h) curve reflects magnetization of the localized- 
moment system. Thus we have a transition from an itinerant (albeit metamagnetic) to 
localized-type behavior clS du function of h. 

Substituting the expressions for m and d 2 taken from Appendix A into the free energy 
(pi]) we can determine the system behavior in the vicinity of h = 0. The phase diagram for 
cases h = and 0.01 , for different values of J/U is exhibited in Fig. 7. Upon increasing of 
J/U the paramagnetic insulating (PI) phase expands at the expense of the paramagnetic 
metallic (PM) phase . In the applied field the boundaries shift towards lower temperatures. 
We see also a typical reentrant metallic behavior at high temperatures. Namely, with rising 
temperature the system evolves from a metal through an insulator back to the metallic 
state. In the inset we display the temperature dependence of the free energies for PM and 
PI states. 

The shape of the phase boundary is essentially the same as in nondegenerate-band case 
||. However, it is shifted remarkably towards lower values of U/W already for rather small 
values of J/W. The boundary is of first order apart from the point specified. The upper 
part of the curve is only qualitative, particularly if realistic DOS is used, as higher-order 
contribution in r will become important. In fact, the solution of the Hubbard model in 
infinite-dimension limit |T(J provides only a crossover behavior, not a weakly discontinuous 
retrograde behavior. 

One may ask if the low-T analysis is realistic, since in the mean-field slave-boson analysis 
the low-energy spin fluctuations (quantum Gaussian fluctuations around the saddle-point 
here) have been neglected. Those spin fluctuations lead to the contribution ~ T 4 In (T/6) 



12 



in the free energy fL4j| . This contribution is of higher order than the T 2 contribution com- 
ing from the quasiparticle excitations across the Fermi surface. Therefore, the analysis is 
realistic, but only to the leading order in low-T expansion, as we have done it. 



B. Physical discussion 

The evolution of the almost localized fermions discussed above for the half-filled band 
case can be explained nicely from a physical point of view. Namely, for U — > U c the 
renormalized band ($e) and the correlation (Ud 2 ) energies almost compensate each other. 
In such situation, much smaller entropy (~ TS) or applied field (~ HbHq) energies tip the 
balance towards either M or I phase. Explicitly, at low temperatures, the spin-disordered 
magnetic insulator has much larger entropy contribution (— fc^Tln2, per orbital) than that of 
the almost localized Fermi liquid (— 7 T 2 / (2$) ,with 7 = (2/3) ir 2 k\p). This circumstance 
tips the balance from PM state (with E@ < 0, but small entropy contribution) towards PI 
state (with Eg = 0, but much larger entropy contribution). At much higher temperatures the 
balance is tipped back towards PM phase, since eventually the entropy of the metallic state 
grows and approaches the asymptotic value 2k b hi 2 per orbital. Thus the reentrant metallic 
behavior is driven by the entropy. It is observed in both (Vi-xCr x ) 2 3 [JT3] and NiS 2 -xSe,j; 
systems. It can be applied also to explain the low-temperature reentrant liquid behavior 
in liquid 3 He ||16|| . In our view the reentrant behavior appearing either as a crossover or as 



a discontinuous transition ( 3 He) uniquely present in the Mott-Hubbard systems defined as 
systems, for which the band and the Coulomb parts of their energy almost compensate each 
other. The exchange contribution tips the balance further towards the localized state. 

Obviously, the first-order nature of the transition will lead to the coexistence of the two 
phases, with localized and itinerant electrons, respectively. However, this mixed phase can 
be discussed only when the magnetism is included, and we will not elaborate on it here. 

As mentioned at the beginning, the present formulation represents for H a = T = a 
simplified version of the full slave-boson and Gutzwiller treatments |IJ discussed recently. 
In Fig. 8 we have compared the Hartree-Fock (HF), ours (SB-HF) and full slave-boson (SB) 
results for the ground state energies for h = and n = 2. the arrows denote the position of the 
Mott-Hubbard boundary in the two latter schemes. The energy difference diminishes rapidly 
with growing h. The difference for physical quantities is only quantitative, not qualitative, 
though the SB approach (]]] has slightly lower energy, since it contains many more variational 
parameters. In the next Section we extend the main features of our solution to arbitrary 
J/W and higher degeneracy d of equivalent orbitals and an arbitrary filling. 

IV. METAL-INSULATOR TRANSITION FOR ARBITRARY ORBITAL 
DEGENERACY AND FILLING: A SPIN ROTATION INVARIANT MODEL 

A. Global (site) representation of the intraatomic interaction. 

We now generalize the principal features of our argument to the case of orbital degeneracy. 
First, we represent the intraatomic part in terms of global (site) representation. For that 
purpose we start from the following expression of that part 
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^ ill aa ill 

where now I and I assume the values 1, 2, . . . , d , and the primed summation is taken for 
I ^ I . We introduce the global spin and particle number operators through the relations 

d 

s t = {s+,s-,s*) = J2 s ^ ( 6 °) 

1=1 

and 

m = ^ n ™ = ^2 nila • ( 61 ) 

a la 

We have the following relations between global and previously introduced operators 

n a n u , = n 2 - m - 2 ^ ^t^U , (62) 
u' i 



and 



^2 'Su ■ s a > =st + ^Yl n m n m ~ \ n * ■ ( 63 ) 
u' i 



Hence, up to constant Hi takes the form 



H I = ^Kj2n 2 i -jJ2s 2 i+lJ2 n ^ n ^ ( 64 ) 

i i il 



3 
2' 

For particular example of e g band we have 



where I — U ' — K — | J. 



^i = l(u-lAj2 n "- J J2 S " + J Yl WA- (65) 

^ " i i il 

The first term represents charge fluctuations, the second the atomic Hund's rule. More- 
over, unlike in a nondegenerate band the intraorbital intraatomic (Hubbard) parameter (J) 
is much smaller than U. Additionally, the first two terms are proportional to d 2 while the 
third is ~ d. Hence, the last term can become the smallest, particularly for highly degen- 
erate systems. In general, the Hamiltonian (64) expresses the so-called minimum polarity 
model of Van Vleck (charge fluctuations suppressed by growing U), as well as the separa- 
tion of the dynamic processes into inequivalent charge and spin degrees of freedom in the 
effective one-band model introduced by Hubbard |I7 |. For highly degenerate systems we 
can approximate 

il il i 
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Noting that 



n i\ n a = -7- G"* • s i) > ( 67 ) 



where /ij is unit vector along an arbitrarily oriented spin quantization axis for the spin Sj, 
one has finally 

^4S^- J S S «^E^-^. (68) 

iii 

where U = 2K + I/d. For a particular case of a nondegenerate band (d = 1), K = J = we 



recover the earlier results [[Lq] with U = U, and I = U. For d 00 the last term is absent 



and the correlated paramagnetic state at T = is described by two variational parameters: 

A = (n?>, (69) 

and 

m={S>). (70) 

Note that in this Section m represent the local-moment magnitude. The third parameter 
describing the long-range order is obtained by making the Hartree-Fock approximation 

• Si) 2 ^ ^{fit- (Si)) fH- S t . (71) 

i i 

The collinear magnetic ordering is expressed then through (Si?) = fii ■ (Si). Since in this 
Section we consider only a paramagnetic state at T = H a = 0, we neglect in our analysis 
the last term. Note however, that the present parameters A, m and (fii ■ Si) for arbitrary d 
correspond directly to the parameters d 2 , pi and p 2 in the preceding Section (e 2 is removed 
via the completeness condition e 2 + p 2 + p 2 + d 2 = 1). 



B. Magnitude of local-moment and charge fluctuations: first-order transition to the 

insulating state at T = 0. 

In direct analogy to the doubly degenerate case discussed in Section II, we express the 
ground state energy in the form 

^ = $(A,m)6 + ^A-Jm, (72) 

where now the band narrowing factor $ depends on A and e represent the average band 
energy for a degenerate system. Without a loss of generality one can assume that $(A, m) = 
A (m) G (m) . In accordance with our simple derivation of the mean-field (Gutzwiller) ap- 
proach in the half-filled band case []T!| we make an expansion: 



A (A) = l + h\ + l 2 X 2 , (73) 



15 



and 



G (m) = g + gim + # 2 ra 2 



(74) 



Note that ^ $ = AG ^ 1. The expansion has the meaning of a Landau expansion, 
and the coefficients can be determined by calculating E G explicitly in limiting situations. 
For example, in the Hartree-Fock approximation we have A — G — l and thus elementary 
analysis provides us with 



A = A = n + n 2 



1 

1 ~ 2d 



and 



n 



(75) 



(76) 



where n is the band filling (n = d corresponds to the half-filling). Analogously, for U — > oo 
and J — > oo, the quantities An and m n will approach their atomic values A™, m m 



oo, the quantities An and m will approach their atomic values Aqo, 

|n for n ^ 1 
f (f + 1) forn^l, 



m = moo 



and 



A = A, 



n for n ^ 1 
n 2 for n > 1. 



Additionally, applying the equilibrium conditions 

dE G dE G 



0, 



(77) 



(78) 



(79) 



we obtain the expression 



A = A 



m = m 



U 



1 + 



8/ 2 A |e|G(m) 
J 

2g 2 m |e| A (A) 



(80) 
(81) 



As in nondegenerate case, the charge fluctuations are suppressed with increasing U. 
On the contrary, the magnetic moment grows with increasing J. Additionally, using the 
conditions A (An) = G (m ) = 1 we obtain that 



A (A) =1 + / 2 (A-A ) 2 , 
G (m) = 1 + g 2 (m - m y 



(82) 
(83) 



The coefficients l 2 and g 2 are determined from the condition that for any integer band 
filling n > 1 we have that A (Aoo) = 0, and G (m^) = n a /n = 1/2 (this result is valid for 
the paramagnetic configuration only and reflects the Pauli exclusion). So, finally we have 
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A(A) = 1 



and 



G(m) 



1 

1 - - 

2 



A-Ap 
Aqo — Ac 



m — mo 



(84) 



(85) 



Substituting these expressions to the expression ( f72[ ) for Eg, we obtain the explicit form 
in the terms of variational variables A and m. Making use of the conditions dEc/dm = 
dEc/dX = we arrive at the algebraic equations, which may be transformed into the 
following equations for the functions G and A taken at the extremal points: 



G s -G 2 



U 



1 

92 



and 



A - A 



l+2 { + (w ) { 



+ A- — 




(86) 



Jr. V J t 



(87) 



with Uq = 8 |e| , and J c = 2 |e| . These two equations can be transformed into each other by 
changes A <-> G, g 2 <-> ^2, and J/J c <-> —U/U . Hence, it is sufficient to solve numerically 
one of them and adapt it subsequently for the second equation. One should notice that the 
nondegenerate-band-case result A = Aoo (i.e. d 2 = 0) is recovered for U/U = 1. 

Before going into the numerical analysis let us summarize the above Subsection. We 
have developed a relatively simple scheme of calculating the magnitudes of local moment 
and of the charge fluctuations (n 2 ), which is equivalent to Gutzwiller-Brinkman-Rice 
scheme || for d = 1. These magnitudes are calculated from a relative balance between the 
renormalized band energy from one side, and the correlation energies ( jA 2 — Jm 2 ) from 
the other. The method involves an interpolation between low- and high-correlation regimes 
(note that the mean-field slave boson theory requires such an interpolation to the U — > 

limit mm). 



C. Numerical analysis for arbitrary degeneracy and filling 

We define the reduced variables 

A « = r^- < 88 > 

and 

m R = °-. 89 

nioo - m 

In this manner, the values \r = = corresponds to the Hartree-Fock approximation, 
whereas the limit with Aj? = m# = 1 correspond to the exact atomic limit, which for an 
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integer n ^ 1 corresponds to the Mott-Hubbard insulator. In Fig. 9 we have plotted Eg as 
a function of Ur = U/U = U/8 |e| , for Jr/Ur = 0.1, which corresponds to J/C7 = 0.4. The 
bare band energy for the featureless density of states is e = — {W/2) n (1 — n/2d) . One sees 
that .Eg = for a critical value of C/r. At this point the values Aoo and are reached 
for n > 1 in a discontinuous way, as illustrated in Figs.lOab. The localization threshold 
diminishes with n. This is because the band energy varies roughly ~ n, while the interaction 
energy is ~ n 2 . The transition to the localized- moment state is determined by the interplay 
between the exchange and the Coulomb interaction. This is illustrated in Figs. 1 lab, where 
the band-narrowing factors A and G have been specified. In most situations the part G (m) 
changes from its Hartree-Fock value only a little. This provides an a posteriori justification 
of the Hartree-Fock approximation for the exchange term. The dominant role of the term 
~ n 2 over the Hund's rule term gives some support to the interpretation of the local-moment 
formation in the terms of a nondegenerate Hubbard model. This becomes clear if one notices 
that in that case 

A = (n 2 ) =n + 2(n^n iL ) = n + 2d 2 , (90) 

and hence A = = n correspond to the limit d 2 = (A = 0). In the situation depicted in 
Fig. 11a A ph 0.6 at the transition. However, unlike in nondegenerate system, the spin-charge 
fluctuation coupling is leading to the discontinuous character of this transition for n > 1. 

The discontinuous nature of the metal-insulator transition induces only a weak enhance- 
ment of the effective mass close to the transition (the curves in Fig. 10b are physically mean- 
ingful only below the discontinuity points). The results in this respect are valid universally 
for arbitrary d > 1 and n > 1. 



V. CONCLUSIONS 

In this paper we have put an emphasis on the similarities and differences of the Mott 
localization in degenerate-band system with the extensive analysis of MIT existing for the 
nondegenerate band case. For that purpose we have made in the first part the Hartree-Fock 
approximation for the exchange interaction. This approximation is applicable in the limit 
when J is substantially smaller than U, as is usually the case for 3d bands. This scheme 
provides us with the physically plausible conclusion that the metal-insulator transition is 
mainly driven by the intraatomic Coulomb interaction; the intraatomic exchange is respon- 
sible for the first-order nature of the transition already at T = 0. Also, with this approach 
one can see that the correlated systems, whether orbitally degenerate or not, can in the 
mean-field approximation be described as systems of fermionic quasiparticles with spin-split 
masses (in the magnetically polarized state) and a nonlinear molecular field coming from 
the correlations, in addition to the usual exchange field. 

The fundamental question is whether the present approach (as well as those listed in 
Ref.l) provide a proper mean- field theory of a correlated state near the metal-insulator 
transition. It seems so and the proper order parameter in the orbitally nondegenerate para- 
magnetic system is either <3> or d 2 , which are nonzero in the metallic phase and vanish in 
the insulating state. What is more important, since for n = 2 the band-narrowing factor $ 
in PM state can be directly related ptl to the physical quantity Z{= $) representing the 
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discontinuity of the Fermi-Dirac distribution at the Fermi energy, the order parameter is a 
measurable quantity for this quantum phase transition at T = 0. In the magnetically ordered 
state the system is additionally characterized by a staggered moment (S?). In the second 
part of this paper we have shown that even in the paramagnetic state the orbitally degener- 
ate system is additionally characterized by the parameter (Sf) describing the local moment 
magnitude. Thus, in our view, the set of the parameters: d 2 <=> (nf), (Sf), and /ij • (Si) 
compose a minimal set describing the metal-insulator transition in the degenerate system 
and associated with it magnetic transition (in nondegenerate system (Sf) = | (1 — 2d 2 )). 
Obviously, the full slave boson and Gutzwiller approaches provide an essentially the 
same qualitative picture, although contain more parameters, which are eliminated by imple- 
menting the constraints appearing as consistency conditions. In this respect, our simplified 
approach provides a didactical guidance for more complicated analysis. For example, the ef- 
fect of quantum Gaussian fluctuations in auxiliary Bose fields neglected so far will introduce 
intersite exchange interactions, which will not introduce any additional order parameter, 
though the detailed thermodynamic properties will contain the contribution coming from 
the interaction between the quasiparticles with the novel characteristics (q a , j3, j3 ), which 
appeared on the mean-field level. 

The slave boson approach has been recently extended [2~2| to describe the antiferromag- 



netic phase in the doubly degenerate Hubbard model. However, the temperature dependence 
of MIT in the paramagnetic phase is still important, since in NiS2_ a: Se a ; || one observes a 
transition to the semiconducting phase upon heating the system, when system crosses the 
Neel point (cf. Fig. 7 for H a = 0, where M— >I transition is observed upon heating the 
system!). The detailed analysis requires the discussion of antiferromagnetic insulating and 
metallic states at T > before any direct comparison with the experiment is made. 
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APPENDIX: MAGNETIC SUSCEPTIBILITY AND DOUBLY 
OCCUPANCY OF THE METALLIC STATE CLOSE TO THE 
MOTT-HUBBARD LOCALIZATION 



We start from the expression (p4]) for the free-energy functional for a constant DOS and 
in the half-band case: 

f = — — = - (4 - m 2 ) $/8 + 2nd 2 - hm - -jm 2 - — r 2 + fc. (A.l) 



where $ is defined by eq. (p4|). 

Minimizing above expression with respect m we obtain: 

df d 2 m , 1 2tt 2 

dm ^4(1 -2d 2 ) 2 -m* 4 3$ 

where $ m = <9$/<9m. Differentiating the above eqation with respect h and taking h = 0, for 
which m = 0, we have 

d 2 I , 8n 2 r 2 _„ 47T 2 r 2 _ 2I . , 

1 _ 2rf2 ^o - 2 - ^0 — ®m$h® X + — g— $mfc$ 2 | = 0, (A.3) 

where m = (dm/dh) h=0 The derivative $ for m = /i = is: 

<9$ f 32d 2 A „ 1 
9m 



(4-m 2 ) 2 V 2 
8rf 



1 - 2d 2 + - \/4 (1 - 2d 2 ) 2 - m?j m + 

=m > (A.4) 



(4-m 2 )v/4(l-2d 2 ) 2 -m 2 



h=0 



We se that $ m= o = 0. Also, (d 2 ) = 0; this is because both $ and d 2 depends only on 
even powers of either m or h. In that situation: 



<9 2 $ 



dmdh 



o 



2d 2 (l - 2c/ 2 + Vl - 4c/ 2 + 4c/ 4 ) m' c 



d 2 



n m n 

VI - 4c/ 2 + 4c/ 4 ° 



4c/ 2 (l-2c/ 2 ) 2 -c/ 2 , d 2 (Ad 2 - 3) (Ad 2 - 1 



m = — m (A.5) 



1 - 2d? u 1 - 2c/ 2 

So, we obtain the magnetic susceptibility x = dm/dh in the form: 



c/ 2 1 , 7T 2 (4ci 2 - 3) (4c/ 2 - 1) 
I -2d 2 ~ 2 J + 48c/ 2 (1 - 2d 2 ) 3 



m o ~ 2 ( - — r — -j + ^— — ^—-3 V 2 . (A.6) 



We use also the low-temperature expansion of ci 2 derived earlier 

1 / u\ T 2 7T 2 U 



d 2 = -[!--) : " - (A.7) 

4 V 2^ 6(1 -in 2 ) 2 
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Expanding m' to the first order in r 2 we arrive at the expression for \ i R physical units: 

I 1(5 , , u 2 + Au-A 



X = Xo 



2-u _ ■ 
2+u J 



, 16 2 2 
H TT U 

3 



(4-?i 2 ) 2 (-2 + m + 2j+jm) 2 



(A.8) 



where xo is the susceptibility for noninteracting particles. 

We see that for u > —2 + 2^/2 « 0.83 the susceptibility is always rising with temperature. 
Note that for DOS smoothly varying around the Fermi energy this increase is solely due to 
the correlations and corresponds to the approaching the localization boundary depicted in 
Fig. 7. In the nondegenerate case and for T = 0, the corresponding Brinkman-Rice formula 
reads: 



X = Xo: 



2 + u 



u 



(A.9) 



Now we derive d 2 in small magnetic field and temperature. Because of time reversal 
symmetry d 2 depends only on even power of h. We determine nonvanishing term ~ h 2 at 
T = 0. Substituting into eq. (fi~2"l) first the expression for m in the first order of h given by (|44"D 
and subsequently differentiating twice, thus obtained formula we see that doubly-occupancy 
probability diminishes with growing T and/or h according to: 



d 2 = 1 -(l-^ 
4 V 2 



7l 2 U 



6(i- tr 



(-2+u+2j+juY 



r/r 



(A.10) 



The decrease of r\ with T and h signals a precursory localization effects discussed in 
subsection II. E. 

We now generalize the results to the case with arbitrary DOS. We start from the expres- 
sion (51) for the free energy with the condition 



E 



W/2 



W/2 



p (e) n a (e) de, 



(A.ll) 



for the chemical potential; n a (e) is given by (|22"D . Defining H = h + jJm — [3, we obtain 
the following expression of m to the first order: 



m 



W/2 



( 



W/2 



Me 



1 + e k B T 



+ 



1 



\ 



l +e k B T 



k B T 



■H 



de 



(A.12) 



or, explicitly: 



with 



m = 2(h + -Jm- f3)K, 



K = 



W/2 



e k b j 



W/2 1 + e 



- k B T 



p (e) de . 



(A.13) 



(A. 14) 
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Differentiating the above equation with respect h and taking h = we have: 



2[l + -Jm -f3 h=0 



K(h = 0) 



Differentiating the free energy with respect to m we obtain the relation: 



and hence: 



i ^ r w/2 ( i \ 

- Jm - (3 + 2 J P ( e ) - a- J) n a (e) = 0, 

rW/2 

/3 = 2$ m V / ep{e)n a (e), 

Z J-W/2 



W/2 

which, when differentiated with respect to h, leads to: 

-W/2 



/W /A 
ep (e) n (e 
-W/2 



W/2 

where Q m h\o is given by (A. 5). From (A. 15) we have: 



dm 
dh 



K 



h=0 



\KJ + AKe{T) 



rf 2 (4d 2 -3)(4rf 2 -l) 



+ 1 



(A.15) 



(A.16) 



(A.17) 



(A.18) 



(A.19) 



2"" 1 V ' l-2d 2 

where all quantities on right hand side are taken for h = 0, and e (T) is average band energy 
per site and orbital. 

Next, we expand e (T) , d 2 and K in the powers of T. From [2U| we have that: 



d 2 



1 2ti 2 Ip/k B T 

4 1 J 3[/„l<P„ 



$ = $ 



o 



and 



where / = U/Uo, and Uo = 8 |e| . Defining //* = yu/$, and T* = T/$, we can write 

^w/2 9 



(A.22) 
(A.23) 

(A.24) 



A" 



1 + e ^ 



■p (e) o?e 



2 



w/2 

/' 

■W/2 

W/2— jt* 



P (e) rfe 



$ I -W/2- 



f (x) p (xk B T* + /i*) c£r 



$2^w J (2n)! 



x 2n f (x)dx 



(A.25) 
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where /' (x) is the derivative of the Fermi-Dirac function. Noting that 

" df (x) 



00 
+00 

x 



dx 
2 df(x 



-dx = — 1, 
1 



dx = --7T 2 , 

dx 3 



we arrive at the expression 

K = 

From |20H we also have: 



2 



pW + ^p (p*)(k B T*f 
o 



^ (T) = e F -^ r (^n 2 ^ ) 
o p 



(A.26) 



(A.30) 



(A.31) 



and thus 



K 



with r 



,(1) 
p 



2p 
$ 

2p 

G (2) 



7T 

~ "6 

+ ■ 



k B T 

167T 2 I 2 

3$g C/ 



(i) 



p(k B T 



P 

*\2 



,(2) 



7T 



6<Pn 



(A.32) 



. We substitute (A.22), (A.24) and (A.32) into (A. 19). In result, the 



susceptibility x — ^Xo {dm/dh) h=0 is of the form: 



X 



Xo 
S 



1 — a 



it 2 (k B Tf 
6$ 2 S 



where 



and 



S = $r 



Up 



1 + 1/2 
(1 + /) 2 



u 2 (i 2 + 21 + 3) (i + iy 2 p 2 - 



jp, 



m 2 P 

Uc 



(A.33) 



(A.34) 



(A.35) 



Neglecting higher-order contribution in J, the result reduces to the usual Stoner form 
X (T = 0) = Xo [1 — P (U + J)]~ . Formula (A.33) generalizes the result of Brinkman and 
Rice H obtained for a nondegenerate case at T = 0, and its generalization [2G] to T > 0. The 
Hund's rule coupling enhances the susceptibility and can be the source of a ferromagnetic 
instability. 

With the help of the expression for $ one can determine the first nontrivial contribution 
to the specific heat in the applied field. Namely, one has C v = (70/^0) T + bH 2 T, where b is 
a constant divergent at the localization point. This increase is connected with the narrowing 
of the band with the growing field ($ ~ d 2 ). 
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Figure Captions 



Fig.l Field dependences of double occupancy d 2 per orbital (a) and magnetization m per 
atom (b), both for T = and n — 2. 

Fig. 2 Jump of the magnetization Am per atom and the critical magnetic field h c both as 
a function of exchange interaction J/W and calculated at the metal-insulator transition for 
T = and n = 2. 

Fig. 3 Regimes of U/W and J/W for the half-filled band configuration comprising both first- 
and second-order transitions regimes in an applied magnetic field. 

Fig. 4 Field dependences of magnetization m per atom, double occupancy d 2 per orbital, and 
the spin-dependent mass enhancement m a /m , all for n < 2 (from the top to the bottom, 
respectively). 

Fig. 5 Regimes of U/W and J/W of first- and second-order transitions in applied magnetic 
field for T = 0. 

Fig. 6 Field dependences of the magnetization and double occupancy for different values of 
J/U. 

Fig. 7 The paramagnetic metal (PM) - paramagnetic insulator (PI) phase boundaries as a 
function of U/W for n = 2 (thick lines), and their shift in the applied magnetic field (fine 
lines). The inset displays the free energies for PM and PI states as a function of temperature. 

Fig. 8 Comparison ground state energies for T = h = 0, n = 2 and J/U = 0.1 obtained from 
the methods: Hartree-Fock (HF), slave boson(SB) and slave boson combined with Hartree- 
Fock (SB-HF). The arrows indicate the position of MIT in the two last approximation 
schemes. 

Fig. 9 The ground energy per site (in units of |e|) as a function of the Coulomb interaction 
Ur = U/Uq = U/8 |e| . A transition to the Mott-Hubbard insulating state takes place for 
each (integer) band filling n for the value of Ur at which Eq = 0. 

Fig. 10 a) The relative magnitude of the charge fluctuations Xr — (A — Ao) / (Aoo — Ao) 
versus Ur. Note a discontinuous nature of the transition to the atomic configuration A = A^ 
at a critical value of Ur specified for each n > 1. 

b) The relative magnitude of the local moment itlr = (m — itlq) / (moo — mo) , a s a 
function of Ur. The discontinuous jump in m reflects the same behavior as of Xr. 

Fig. 11 a) The band-narrowing parts A (A) and G (m) as the magnitude of the Coulomb 
interaction. The inset compares the two narrowing factors, both diminishing with increasing 

Ur. 
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b) The charge- and spin-fluctuation parts of the band narrowing as a function of 
the Coulomb interaction, for the value of Ur below the Mott transition. The inset displays 
the different trend as a function J /If ratio. 
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